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Abstract 



We improve the calculations of the n — ► 7r°77 decay within the context of me- 
son chiral lagrangians. We use a chiral unitary approach for the meson-meson 
CNJ . interaction, thus generating the ao(980) resonance and fixing the longstanding sign 

ambiguity on its contribution. This also allows us to calculate the loops with one 
vector meson exchange, thus removing a former source of uncertainty. In addition 
CN . we ensure the consistency of the approach with other processes. First, by using 

vector meson dominance couplings normalized to agree with radiative vector meson 
■ decays. And, second, by checking the consistency of the calculations with the re- 

lated 77 — > tt°t] reaction. We find an n — > 7r°77 decay width of 0.47 db 0.10 eV, in 
clear disagreement with published data but in remarkable agreement with the most 
recent measurement. 



§^ 1 Introduction 



The r\ — > 7r°77 decay has attracted much theoretical attention, since Chiral Perturbation 
Theory (ChPT) calculations have sizable uncertainties and produce systematically rates 



about a factor of two smaller than experiment jU EJ. In contrast models using quark 
box diagrams [3J E] claim to obtain acceptable rates. Within ChPT, the problem stems 
from the fact that the tree level amplitudes, both at 0(p 2 ) and 0(p 4 ), vanish. The first 
no n- vanishing contribution comes at 0(p A ), but either from loops involving kaons, largely 
suppressed due to the kaon masses, or from pion loops, again suppressed since they violate 
G parity and are thus proportional to m„ — jHj. The first sizable contribution comes 
at 0(p 6 ) but the coefficients involved are not precisely determined. One must recur to 
models: either Vector Meson Dominance (VMD) 00171, the Nambu-Jona-Lasinio model 
(NJL) jH], or the extended Nambu-Jona-Lasinio model (ENJL) [3110], have been used 
to determine these parameters. However, the use of tree level VMD to obtain the 0(p 6 ) 
chiral coefficients by expanding the vector meson propagators, leads |5J to results about 
a factor of two smaller than the "all order" VMD term, which means keeping the full 
vector meson propagator. All this said, the chiral approach has been useful to unveil the 
physical mechanisms responsible for this decay, but it has become clear that the strict 
chiral counting has to be abandoned since the 0(p 6 ) and higher orders involved in the full 
("all order") VMD results are larger than those of 0(p 4 ). For a review, see jTT], together 
with an experimental upper bound. 
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Once the "all order" VMD results is accepted as the dominant mechanism, one cannot 
forget the tree level exchange of other resonances around the 1 GeV region. In comparison 
with VMD, the exchange of J PC = l + ~ axial vectors f21EI] yields negligible contributions 
when using values of the couplings in agreement with 77 — > 7r°7r° data ^3]. Still at tree 
level, the ao(980) exchange, which was taken into account approximately in j^j, was one 
of the main sources of uncertainty, since even the sign of its contribution was unknown. 

After the tree level light resonance exchange have been taken into account, we should 
consider loop diagrams, since meson-meson interaction or rescattering can be rather 
strong. First of all we find the already commented 0(p 4 ) kaon loops from ChPT, but 
also the meson loops from the terms involving the exchange of one resonance. The uncer- 
tainty from the latter was roughly expected [5] to be about 30% of the full width. 

Another relevant question is that no attempts have been done to check the consistency 
of n — > 7r°77 results with data from other processes. On the one hand, the decay results 
have not been compared with the crossed channel 77 — > 7i°r), although some consistency 
tests with 77 — > 7T 7r° have been carried out as quoted above. The reason is not surpris- 
ing since there are no hopes within ChPT to reach the ao(980) region where there are 
measurements of the 77 — > tt i] cross section ^SJ E>] . On the other hand, the explicit 
SU(3) breaking already present in the radiative vector meson decays has not been taken 
into account when calculating the VMD tree level contributions. 

The former discussion has set the stage of the problem and the remaining uncertainties 
that allow for further improvement. In recent years, with the advent of unitarization 
methods, it has been possible to extend the results of ChPT to higher energies where the 
perturbative expansion breaks down and to generate resonances up to 1.2 GeV [T71 FTSl 
HH 1201 El I22j • In particular these ideas were used to describe the 77 — > meson — meson 
reaction, with good results in all the channels up to energies of around 1.2 GeV [231- Work 
in a similar direction for this latter reaction has also been done in |2*H |23 12H1 • With these 
techniques, and always within the context of meson chiral lagrangians, we will address 
three of the problems stated above: First, the a (980) contribution, second, the evaluation 
of meson loops from VMD diagrams and, third, the consistency with the crossed channel 
77 — > 7t°n. In particular, we will make use of the results in [23], where the 77 — > 7T°n 
cross section around the ao(980) resonance was well reproduced using the same input as 
in meson meson scattering, without introducing any extra parameters. 

With these improvements we are then left with a model that includes the "all order" 
VMD and resummed chiral loops. We expect this approach to provide a good description 
of 77 — > 7r°77 because recent studies on the vector meson decay into two pseudoscalar 
mesons and one photon fZT\ I2H1 120] indicate that such a combination of " all order" VMD 
contribution plus the unitary summation of the chiral loops leads to good agreement with 
data in a variety of reactions. These include <fi —* 7r°7r°7 |29j . where the chiral loops are 
dominant, uo — > 7r°7r°7 [23I2H], where the VMD mechanism is dominant, and p — > 7r°7r°7 
[23 |2H1 where both mechanisms have about the same strength. 

Concerning the fourth issue of the SU(3) breaking present in radiative vector me- 
son decays, we will take it into account here by using effective couplings normalized to 
reproduce the most recent experimental data. 

Incidentally, there are preliminary results from a very recent experiment [20]; which 
give a decay width about half of the previous one. In that work the authors refined 
the background subtraction, which was known to be rather problematic. Let us remark 
that, in view of the former discussion, revisiting the previous theoretical calculations is 
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mandatory regardless of whether these new experimental results are confirmed or not. 

In what follows we will address all these theoretical issues in detail, including an 
updated estimation of the uncertainties in the calculation. In particular, we will take 
into account the experimental errors in the radiative vector meson decays, which were 
neglected before, although they will turn out to be the largest source of uncertainty. 



2 VMD contribution 

Following [5 we consider the VMD mechanism of Fig. 1 which can be easily derived from 
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Figure 1: Diagrams for the VMD mechanism. 

the VMD Lagrangians involving VVP and V7 couplings 
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where and P are standard SU (3) matrices constructed with the nonet of vector mesons 
containing the p, and the nonet of pseudoscalar mesons containing the n, respectively. 
For instance for pseudoscalar mesons P = P + with P given by [SI] 
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and similarly for vector mesons jozj. We also assume the ordinary mixing for the (f), the 
u, the 77 and rf: 
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and / = 93 MeV, with Gy the coupling of p to 
From Eq. ((TJ) one can obtain the radiative widths for 



k 3 



(4) 



where k is the photon momentum for the vector meson at rest and Ci are SU(3) coefficients 
that we give in Table 1 for the different radiative decays, together with the theoretical 
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Table 1: SU(3) coefficients together with theoretical and experimental branching ratios 
for different vector meson decay processes. 



(using Gy = 69 MeV and / = 93 MeV) and experimental [2] branching ratios. We shall 
refer to these results as those with "universal couplings". 

The agreement with the data is fair but the results can be improved by incorporat- 
ing SU(3) breaking mechanisms |34j . For that purpose, we will normalize here the C, 
couplings so that the branching ratios in Table 1 agree with experiment. These will 
be called results with "normalized couplings". In this way we are taking into account 
phenomenologically the corrections to the VP'j vertex from an underlying field theory. 

Once the VP7 couplings have been fixed, we can use them in the VMD amplitude 
corresponding to the diagram of Fig. 1, which is given by 
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where q = P — k±, q' = P — k 2 , with P, ki, k 2 the momentum of the 77 and the two photons. 
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We have parametrized the p width phenomenologically as: T p (q, s, s __ 
whereas for the u we have considered a constant = 8.44 MeV. Nevertheless, our results 
are rather insensitive to these details. From the above amplitude, the n decay width is 
easily calculated, as well as the 77 invariant mass distribution 
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where we take for reference the momentum of the pion, p, in the z direction, so that 
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with ui the energy of the ir°. 
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In Fig. 2 we show the results of the mass distribution with and without the radiative 
widths normalization factors. The integrated width is given by T = 0.57 eV (universal 
couplings); T = 0.30±0.06eV (normalized couplings), where the error has been calculated 
from a Monte Carlo Gaussian sampling of the normalization parameters within the errors 
of the experimental branching ratios of Table 1. Let us note that there have been stable 
values for the vector meson radiative widths throughout the last decade in the PDG but 
a sizable change in the PDG2002. Had we used the PDG2000, we would have obtained 
0.21 ± 0.05 eV (normalized couplings). 

It is interesting to compare these results with those in jS], where they used a universal 
SU(3) coupling with Gy adjusted to the uj — > 7r°7 decay data existing at that time, 
and obtained an "all orders" value of 0.31 eV. The difference between that value and the 
0.21 ± 0.05 eV that we would have obtained with older data has to be attributed to the 
adjusting to all the branching ratios, instead of just one as in jS]. 




Figure 2: Invariant mass distribution of the two photons with VMD terms only. The 
solid curve has been calculated with an universal coupling, whereas the dashed one has 
the couplings normalized differently to fit the radiative decays. 

Our VMD normalized result is within three standard deviations from the value presently 
given in P,j2]: T = 0.84±0.18eV, but within one sigma of the more recent one presented 
in [3U|, T = 0.42±0.14eV. There are, however, other contributions that we consider next. 

3 Meson loops 

The contribution of pion loops to r\ — > 7r°77, evaluated in [S], proceeds, to begin with, 
through the G-parity violating i] — > 7r°7r + 7r~ process. Since the contribution is propor- 
tional to m u — ma, it is very small and we think that if such terms are included, other 
isospin violating terms proportional to m u — m^, and isospin violating corrections to the 
main terms should also be included. Rather than undertaking this delicate task, we will 
use the results of j^j to estimate uncertainties from all these sources. 

The main meson loop contribution comes from the charged kaon loops, calculated at 
0(p A ) in El El QUI EE2 d, and proceeds via rj -> tt°K + K- -> vr°77. Note that these 
loops are also suppressed due to the large kaon masses. That is why the r\ — > 7r a (980) — > 



5 



7r°77 mechanism was included explicitly, with uncertainties in the size and sign of the 
a (980) couplings. As commented in the introduction, the chiral unitary approach solves 
this problem by generating dynamically the ao(980) in the K + K~ — > 7T°r] amplitude. 

In this section, we will illustrate this approach by revisiting the work done in [22 on 
the related process 77 — > ir°T] where the chiral unitary approach was successfully applied 
around the ao(980) region. Since for the rj decay the low energy region of 77 — > ir r) is 
also of interest, we will include next the VMD mechanisms also in this reaction. Once we 
check that we describe correctly 77 — > 7r°r), the results can be easily translated to the eta 
decay. We will finally add other anomalous meson loops that are numerically relevant for 
eta decay but not for 77 — > 7i r]. 



3.1 The unitarized 77 —> tt t] amplitude in the ao(980) region 



In [23] it was shown that, within the unitary chiral approach, the 77 — > 7r°?7 amplitude 
around the ao(980) region, diagrammatically represented at one loop in Fig. 3, factorizes 



as 



— it — (t x K + tAK+K-)tK+K-,ir°r] 

with tK+K-,TT°ri the full K + K~ — > TT°r] transition amplitude. 
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Figure 3: Diagrams for the chiral loop contribution 

The first three diagrams correspond to t X KtK+K-,n°ri °f Eq. ©, already evaluated at 
one loop in [313 EH], where the factorization of the leading tx+K-^ri also occurred. In our 
case t X K, written in a general gauge to be also used for the 77 — > 7r°77 reaction, is given 
by 
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above the K + K~ threshold, with the —vk term removed below threshold. Note that 
the unitarized tK+K-,n°r] transition matrix, not just the lowest order chiral amplitude, 
is factorized outside the loop integral. This on shell factorization was shown in (22] by 
proving that the off shell part of the meson-meson amplitude did not contribute to the 
loop integral. 
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The meson meson scattering amplitude was evaluated in [17J by summing the Bethe 
Salpeter (BS) equation with a kernel formed from the lowest order meson chiral La- 
grangian amplitude and regularizing the loop function with a three momentum cut off. 
Subsequently, other approaches like the inverse amplitude method [T[JJ [22] or the N/ D 
method [20J were used and all of them gave the same results in the meson scalar sector. 
For 77 — > 7r°?7 below 1 GeV only the L = 0, I = 1 amplitudes are needed [37]. The BS 
equation sums the diagrammatic series of Fig. 4, which implies that in the 77 — > 7i°r] 



Figure 4: Diagrams summed in the BS equation, using the 0(p 2 ) ChPT vertices, 
transition of Fig. 3 one is resumming the diagrams of Fig. 5. Furthermore, the same 




K 



K 



K 



KK 

7C°T| 



Figure 5: Resummation for 77 — > 7r°r]. 

on-shell factorization of the t matrix in the loops found for 77 — > ir°r] was also justified 
for meson-meson scattering in JJJ. Thus, the BS equation with coupled channels can be 
solved algebraically, leading to the following solution in matrix form 



t(s) = [l-t 2 (s)G(s)]-H 2 (s), 



IF 



with s the invariant mass of the two mesons, t 2 the lowest order chiral amplitude and 
G(s) a diagonal matrix, diag(G-^ x , G V7T ), accounting for the loop functions of two mesons. 
These G functions were regularized in [T7J by means of a cut off. The G analytic expres- 
sions, both using a cut off or dimensional regularization can be found in |2Uj . 

In Eq. flSJ) there is another term, tAK+K-tK+K-,n°r], which corresponds to the last two 
diagrams of Fig. 3 where the axial vector meson i^i(1270) is exchanged. For the one loop 
result we follow |38j . Given the large mass of the axial vector, both the factorization of 
the unitarized on shell meson-meson scattering amplitude outside the loop, as well as that 
of the 77 — ► K + K~ amplitude are also justified [22] ■ Hence, when the full series of Fig. 5 
is considered, one obtains the contribution tAK+K-tK+K-^ ^ with Iak+k- given by 
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First of all we show in Fig. 6 the result for the 77 — > 7i°r] cross section obtained 
from Eq. @, which coincides with that obtained in [23] • To ease the comparison with 
experimental data we also show the events concentrated in bins of 40 MeV, roughly like 
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Figure 6: 77 — > ir r] cross section, using Eq. Z is the maximum value of cos(6) 

integrated. The experimental data come from (SHI EH], the latter ones normalized in 
the a2(1320) peak region. The dashed histogram corresponds to the convolution over an 
experimental resolution of 40 MeV. 

the experimental ones. We can easily notice the peak of the ao(980) whose dynamical 
generation is guaranteed by the resummation of diagrams in Fig. 5. The resummed 
tK+K-^-K°ri amplitude has indeed a pole in the complex plane associated to the a (980) 
resonance |17j . 

In Fig. 6, we also show results above and below ao(980), whose description requires 
further ingredients than those needed just for the ao(980) region. In particular, the 
a2(1320) resonance (second peak) is included phenomenologically as in ref. [23] . 

In [2B! loops like those in Fig. 3, but exchanging a vector meson instead of an axial- 
vector meson were estimated negligible in the do (980) region and hence neglected. In 
addition, the VMD tree level mechanism of Fig. 1 (with an outgoing rj) was neglected 
since it has no resonant structure in the 77 s-channel. As a consequence the agreement 
of Eq. with experiment is fair but some discrepancies can be noticed in Fig. 6 at low 
energies. 

3.2 VMD mechanisms in 77 — ► rjir 

For the purpose of the present work, rj — > 7r°77, the low energy region of the 77 — ► 7i°i] 
reaction is also relevant. Therefore we will include as a novelty both the VMD tree level 
contribution as well as the loops involving vector meson exchange. 

First, we can see in Fig. 7 that the results obtained adding the tree level VMD ampli- 
tude normalized to the u, p radiative decay rates (dashed line) are acceptable around and 
below the peak of the ao(980) resonance. Let us note that the inclusion of these terms 
does improve the description of the low energy region. The binning of the theoretical 
results would make again the apparent agreement with data to look much better, but for 
the sake of clarity we have not added more lines to the figure, as long as the binning effect 
has already been illustrated in Fig. 6. Although in section one we have justified the use 
of the normalized couplings, we also show in Fig. 7 the result using universal couplings 
(dotted line). In this process, the effect of normalizing the couplings of the vector me- 
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son radiative decays is not as drastic as in Fig. 2 for the 77 decay where only the VMD 
mechanism was considered. In what follows we will only use the normalized couplings. 
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Figure 7: 77 — > 7r°i] cross section, using Eq. (0 (continuous line), adding the universal 
VMD contribution (dotted line) or the normalized VMD contribution (dashed line). 



Second, in addition to the axial vector meson exchange in loops considered in the 
previous section, we have to include the loops with vector meson exchange for complete- 
ness. In fact, some of the uncertainties estimated in [5] were linked to these loops. For 
consistency, once again we have to sum the series obtained by iterating the loops in the 
four meson vertex shown in Fig. 8. Hence, the new amplitude, which we shall call t VMDL , 
will be given by: 
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(13) 



where now i^^ 40 is the factor that multiplies the e±e 2 product in the s-wave projection of 



the tL^ D amplitude in the 77 — > 7r°r] CM. Although the Lorentz structure of polarization 



vector products may seem rather complicated from Eq. (0), it is easy to show that after 
the s-wave projection the polarization vectors factorize indeed as t\£ 2 . In a general frame 
the ei€ 2 factor has to be replaced by t\t 2 — (k 2 t\)(k\t 2 ) j {kxk^ . Once again we have 
factorized the amplitudes for the same reasons as done with the other terms. 




Figure 8: Loop diagrams for VMD terms. The diagrams with the two crossed photons 
are not depicted but are also included in the calculations. 
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Of course, when introducing loops with vector meson exchange we have to consider 
loops involving a K* + or a K*° exchanged between the photons (see Fig. 8.b), which were 
not present at tree level. These would be given by 



The contribution of all these new VMD loop diagrams is an increase of the order of 
10-20% of the result shown in Fig. 7 by the dashed line (normalized VMD couplings). 
The new result would overlap in a large region with the dotted line of Fig. 7 and hence 
we do not show it explicitly. 

In what follows we make some considerations about the diagrammatic interpreta- 
tion of the "all order" VMD calculation, the normalization of the VP'y vertices and the 
meson-meson interaction. By "all order" VMD one means [S] that the full vector meson 
propagator (s — My+iM v T(s))~ l is used in the calculations. This full propagator includes 
self-energy diagrams in a Dyson- Schwinger resummation, leading to a shift of the bare 
mass and generating a width j^T]. Thus, one must think in terms of self-energy insertions 
in the middle of the vector meson lines in Fig. 1 and 8. The VP'j coupling normalization 
to agree with the radiative vector meson decays can also be understood as considering 
vertex correction diagrams in Fig.l and 8, and therefore it does not lead to any double 
counting with the dressing of the vector meson propagator. Finally, the meson-meson 
interaction in the VMD terms leads to the diagrams of Fig. 8, in which the two pseu- 
doscalar mesons interact through four-pseudoscalar meson vertices. The resummation of 
pseudoscalar meson-meson loops thus leaves apart the vector meson lines and the VP^f 
vertices. Once again this ensures that there is no double counting. 

3.3 Meson loops in rj — > tt°77 

In the rj —>■ 7r°77 case, the meson loop diagrams correspond to those of n r] — > 77 but 
considering the 7r° as an outgoing particle. Hence, it is enough to replace s = (p v + p n ) 2 
by M 2 = (p v — p-x) 2 = (p 7 i + p l2 ) 2 in all the TT°r] — > 77 amplitudes, which factorize in all 
the loop diagrams that we have considered so far, and in the t x x and t^x+K- function. 

Since we are considering all the VMD diagrams and the chiral loops, we still have to 
take into account another kind of loop diagrams |5], shown in Fig. 9, which involve two 
anomalous 7 — > 3M vertices. Despite being 0(p 8 ) it has been found [S] that they can 
have a non negligible effect on the 77 decay. The further rescattering of the mesons in the 
diagrams of Fig. 9, given the structure of the jMMM vertex [5] in the momenta of the 
particles, would be suppressed by factors of p 2 /q 2 (with q the loop variable) with respect 
to those considered for the VMD mechanism. This, and the fact that these anomalous 
terms are already small, makes the consideration of rescattering in these loops superfluous. 
Therefore, it is enough to take the results from [S] where it is found that their largest 
contribution comes from the kaon loops. We use Eqs.(12),(13) and (27) of that reference 
(note that there is a global change of sign with respect to our notation). Concerning 
77 — > r 7 7r °! these kind of loops have been neglected in the previous section, because the 
intermediate particles are very far off shell, due to the crossed character of the loop in the 
reaction. 
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Figure 9: Diagrams with two anomalous 7 3M vertices. 



4 Results for 77 — > 7r 77 

Using the model described in the previous sections, we plot in Fig. 10 the different contri- 
butions to dT/dMj. We can see that the largest contribution is that of the tree level VMD 
(long dashed line). Let us recall that this is a new result as long as we are using the VMD 
couplings normalized to agree with the vector radiative decays. The resummation of the 
loops in Fig. 3 using Eq.©, (short dashed line) gives a small contribution (0.011 eV in 
the total width), but when added coherently to the tree level VMD, leads to an increase 
of 30 % in the r\ decay rate (dashed-dotted line). More interestingly, the shape of the 77 
invariant mass distribution is appreciably changed with respect to the tree level VMD, 
developing a peak at high invariant masses. The resummed VMD loops in Fig. 8, using 
Eqs.(JT3J) and (HU), leads, through interference, to a moderate increase of the rj decay rate 
(double dashed dotted line), smaller than that of the chiral loops considered before. The 
last ingredient is the contribution of the anomalous mechanisms of Fig. 9 (continuous 
line), leading again to a moderate increase of the 77 decay rate, also smaller than the 
chiral loops from Eq.Q. These anomalous mechanisms have a very similar shape to the 
tree level VMD and interfere with it in the whole range of invariant masses. 
Altogether, when integrating over the invariant mass, we get: 

T(r] -»• 7T°77) = 0.47 ±0.08 eV. (15) 

Note that the inclusion of the loops has increased the tree level VMD result by 50%. 
For comparison, we quote here what we would have obtained using the universal VMD 
couplings: 0.80 eV. 

So far, the theoretical error has been obtained only from the propagation of the exper- 
imental errors in the vector meson radiative decay branching ratios, given in Table 1. The 
errors from this source had not been considered before although they will turn out to pro- 
duce the largest uncertainty. In practice we generate a Gaussian weighted random value 
for each VMD coupling which yield a result for the width. This procedure is repeated 
a sufficiently large number of times, leading to an approximate Gaussian distribution of 
results from where we obtain a central value and the error. 

We come now to revisit the uncertainties considered in [5] . One of the largest sources 
to the ±0.2 eV accepted uncertainties in that work, was the contribution of the ao(980), 
whose sign was unknown. This problem is solved in the present work since the ao(980) is 
generated dynamically from the rescattering of the mesons implied in the Bethe-Salpeter 
resummation of the tx+K-^-K amplitude. Hence its effect can be easily observed compar- 
ing with the standard ChPT result, which is obtained by substituting the full t K + K - vl ,o 
by its lowest order 0(p 2 ). In Fig. 10 this corresponds to the difference between the con- 
tinuous and the dotted line. The contribution of the ao(980) resonance tail is rather small 
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Figure 10: Contributions to the two photon invariant mass distribution. From bottom 
to top, short dashed line: chiral loops from Eq.®; long dashed line: only tree level 
VMD; dashed-dotted line: coherent sum of the previous mechanisms; double dashed- 
dotted line: idem but adding the resummed VMD loops; continuous line: idem but 
adding the anomalous terms of Fig. 10, which is the full model presented in this work 
(we are also showing as a dotted line the full model but substituting the full tx+ k~ ,r)n° 
amplitude by its lowest order). 

and increases the r\ decay rate from 0.47eV to 0.48 eV. The sign of its contribution is 
unambiguously determined. Thus, the present calculation removes completely this source 
of error. The explicit calculation of the ao(980) contribution giving such a small effect 
justifies neglecting the 02 (1320) resonance contribution which lies much further away in 
energy. 

The other source of uncertainty in [Hj was the contribution of the VMD loops. We 
have been able to calculate them in this work and, as seen in Fig. 10, these effects are 
also rather small. They increase the 77 decay rate by 0.02 eV. Altogether the ao(980) plus 
the VMD loops increase the 77 decay rate by 0.03 eV. We thus eliminate these two sources 
of previous uncertainties, while realizing at the same time that the uncertainties of 0.2 eV 
attributed to these sources in jHj were indeed a generous upper bound. 

In our approach the tree level exchange of the /ii(1170), &i(1235) and /ii(1380) axial 
resonances ^2] will be included as an uncertainty. The reason is that, according to |12j . 
they would increase the decay width by about 0.07 eV. However, as shown in f2J|H], their 
inclusion in 77 — > 7r 7r° with the couplings used in ^2] would overestimate the 77 — > 7r°7r° 
cross section. In view of these discrepancies, we thus consider safe to accept a theoretical 
uncertainty of the order of 0.05 eV which should still be a generous upper bound. 

As commented at the beginning of section three, there are uncertainties due to isospin 
violating terms. We estimate the errors from this source using the results obtained in 
|5| for the G-parity violating term corresponding to Fig. 3 but with pion loops. This 
contribution is of the order of 0.05 eV to the total 77 decay rate. 
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Finally, by summing all errors in quadrature, we arrive to 
r(?7 -> 7T°77) = 0.47 ± 0.10 eV 



(16) 



Note that although we have considered a new error source from the uncertainties in 
the vector radiative decays, which turns out to be the largest one, we still have reduced 
the uncertainty from previous calculations. 

The result of Eq. ()16|) is in remarkable agreement with the latest experimental num- 
ber [30 1 . and lie within two sigmas from the earlier ones in [1. 2 . Confirmation of those 
preliminary results would therefore be important to test the consistency of this new ap- 
proach. Furthermore, precise measurements of the 77 invariant mass distributions would 
be of much help given the differences found with and without loop contributions. 

5 Conclusions 

We have reanalyzed the 77 — > 7r°77 decay within the context of meson chiral lagrangians, 
gathering all the mechanisms discussed in the literature, but improving them in the fol- 
lowing aspects: 

On the one hand, using the well tested chiral unitary approach, we have removed the 
uncertainties from the do (980) resonance as well as those from loops with the exchange 
of one vector meson. In particular, since the ao(980) is generated dynamically from the 
meson loop resummation, we have unambiguously determined the sign of its contribution, 
whereas for the one vector loops we have performed an explicit calculation that in previous 
works had only been considered as a large source of uncertainty. 

On the other hand, we have also checked the consistency with other related processes: 
First, a relevant observation is that the tree level vector meson dominance amplitude with 
a universal SU(3) vector-vector-pseudoscalar coupling is not consistent with the p — > 777, 
uj — > 7r°7 and u — > 777 decays. Consequently, throughout the 77 — > 7r°77 calculation, we 
have used couplings normalized to agree with the radiative vector meson decays. Second, 
we have established the consistency of our 77 — > 7r°77 model with the related 77 — > r]n° 
process. 

Furthermore, we have performed a careful error analysis of our results. As a novelty 
we have considered the experimental errors in the vector meson radiative decay widths, 
which turn out to be the largest source of uncertainty. However, since, as just commented 
above, we have removed former sources of uncertainty, our final error is still smaller than 
previous estimates. 

Altogether we have found a result of T(i] — > 7r°77) = 0.47 ± 0.10 eV. 

With the improved calculation just presented, it seems clear that the mechanisms thus 
far suggested in the literature in the context of meson chiral lagrangians lead to a result at 
variance with the experimental result T = 0.84±0.18eV from P3I2]- However, it is worth 
noticing the agreement of the above result with the new preliminary measurement T = 
0.42 ± 0.14eV from [30J. Nevertheless, a measurement of the invariant mass distribution 
would be more stringent. Confirmation of the preliminary results of [30 and an accurate 
measurement of the 77 invariant mass distribution should be the experimental priorities 
to clarify the situation. 
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